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I.  Introduction 

Solvation  of  aromatic  molecules  by  simple  hydrocarbons  can  be  readily 
explored  by  supersonic  molecular  jet  spectroscopy.  The  study  reported  in 
this  paper  expands  the  range  of  systems  investigated  to  N-heterocyclic 
aromatic  molecules. 

We  have  recently  reported  several  studies  of  van  der  Waals  (vdW)  clusters 
in  which  cluster  geometry,  binding  energies  and  the  cluster  nucleation 
processes  have  been  elucidated. ^ The  data  in  these  instances  were 
obtained  using  fluorescence  excitation  (FE),  dispersed  emission  (DE),  and 
two-color  time  of  flight  mass  spectroscopy  (2-color  TOFMS).  Coupled  with  the 
spectra,  potential  energy  calculations  using  atom-atom  potentials  of  the 
exponential-six  (exp-6)  or  Lennard-Jones  (LJ)  form  can  provide  further 

7 

detailed  insight  into  cluster  geometries  and  binding  energies.  The 
calculations  provide  information  that  complements  the  spectroscopic  data  and 


aids  in  the  identification  of  stable  geometries  and  in  the  evaluation  of 
potential  energy  surfaces.  The  constants  employed  in  these  atom-atom 
potential  calculations  are  obtained  through  independent,  published  parameters 
that  are  not  fit  to  our  spectroscopic  results.  The  calculations  thereby  serve 
as  a  discrete  data  set  which  can  be  used  to  verify  the  experimental 
conclusions. 

Our  previous  cluster  studies  have  dealt  with  benzene,  toluene,  aniline. 


and  n-propyl  benzene  solvated  by  the  simple  hydrocarbons  methane,  ethane,  and 
propane.  Significant  progress  has  been  made  toward  a  full  understanding  of 
these  systems^ in  terms  of  binding  energies,  relaxation  dynamics,  and 
geometries.  Two  new  directions  are  now  being  pursued  for  the  extension  of 
these  cluster  studies;  the  study  of  clusters  of  the  above  aromatics  with 


hydrogen  bonding  solvents  (i.e.,  H^O  and  NH^)  and  the  study  of  heterocyclic 
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aromatic  solutes  closely  related  to  biologically  relevant  systems  (i.e.,  DNA 
'  bases  and  amino  acids). 

Pyrazine  (C^H^N2)  is  chosen  as  the  first  solute  system  to  be 

investigated  with  these  future  directions  in  mind.  Pyrazine  excitation 

spectra  in  the  n-»-ir  transition  region  (  ~  3300  X)  are  well  known®  and 

g 

thus  the  study  of  pyrazine  solvation  becomes  readily  accessible. 

Pyrazine  clusters  with  these  hydrocarbon  solvents  are  in  most  respects 
similar  to  the  benzene  hydrocarbon  clusters,  although  subtle  differences  in 
geometry  due  to  the  nitrogen  atoms  are  found.  The  LJ  potential  is  employed 
for  the  calculations  because  many  independent  atom-atom  parameters  are  reported 
for  it;  moreover,  even  specific  hydrogen  bonding  interactions  can  be  incorpor¬ 
ated  into  the  calculations  as  appropriate,  through  the  LJ  potential  in  the  10- 
12  form.''® 

II.  Experimental  Procedures 

The  experimental  procedures  employed  in  this  research  are  similar  to 
those  reported  earlier.'”^  Two  separate  lasers  are  used  to  generate  the 

ions  for  2-color  TOFMS:  the  output  of  the  first  laser  is  scanned  through  the 

# 

n-^  11  pyrazine  absorption  and  the  output  of  the  second  laser  provides  the 
ionization  photon,  taking  the  pyrazine  from  the  ^ state  to  the  positive 
ion,  C^H^N^.  The  pump  laser  consists  of  a  Hd  ^/YAG  laser  pumping  a  DCH 
dye  laser  (  -  648  nm)  whose  output  is  frequency  doubled.  The  ionization  laser 
is  of  the  same  configuration  but  the  dye  laser  operates  with  R610  dye;  the 
frequency  doubled  output  from  the  dye  laser  is  further  mixed  with  1.064  pm 
from  the  Nd^^/YAG  fundamental  (-  227.3  nm). 

The  apparatus  consists  of  two  chambers  separated  by  a  skimmer  and  a  2" 
gate  valve  with  the  (pulsed  or  CW)  nozzle  in  the  first  chamber  and  the  mass 
spectrometer  in  the  second  chamber.  In  this  configuration,  pyrazine  signals 
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are  weak  due  to  the  relatively  low  ( ~  12  torr  @  500  K)  vapor  pressure  of 
pyrazine.  A  new  pulsed  valve,  identical  in  design,  operating  characteristics 
and  operating  conditions,  mounted  in  the  second  chamber,  is  used  for  these 
experiments.  The  mass  spectrometer  ionization  region  and  flight  tube  are 
separated  from  the  nozzle  by  a  skimmer;  the  distance  between  the  nozzle  and  the 
ionization  point  is  about  5  cm.  The  new  configuration  increases  the  signal 
strength  by  almost  a  factor  of  10. 

Special  attention  must  be  given  to  fragmentation  processes  in  these 
experiments.  If  the  ionization  energy  is  too  large,  the  solute-solvent 
clusters  fragment  and  the  mass  selective  information  in  the  experiment  is  lost. 
The  pyrazine  ionization  threshold  is  -44000  cm”^  above  the  first  excited 

singlet  state.  If  as  little  as  1000  cm~^  more  energy  is  employed  to  ionize  the 
clusters,  significant  amounts  of  fragmentation  obscure  and  interfere  with  the 
desired  signal.  Thus,  signals  from  fragmented  C^H^N2(CH^)2  and  C^H^N2(CH^)^ 
can  appear  in  the  C^H^N2(CH^)^  mass  channel  at  intensities  comparable  to  the 
C^H^N2(CH^)^  signal  itself.  Fragmentation  is  not  so  serious  a  problem  for 
the  other  cluster  systems  we  have  studied.^ In  these  studies,  the 
ionization  energy  for  all  clusters  is  kept  well  below  44,000  cm”^ . 

In  conjunction  with  the  experiments,  computer  calculations  are  performed 
using  an  exp-6  or  LJ  potential  energy  function.  These  calculations  give 
stable  cluster  geometries,  binding  energies  and  barriers  to  internal  cluster 
motion.  The  potential  form  employed  for  these  calculations  is 


3 


in  which  the  first  term  constitutes  the  pairwise  electrostatic  contribution 


obtained  from  interaction  between  atom  i  and  j  of  different  molecules, 

and  q.  are  partial  electronic  charges  (in  electronic  units)  centered  on  each 

5  •  -1  / 

atom,  and  1.16  x  10  is  a  conversion  factor  to  give  E  in  cm  if  r.  .  (the 

^  J 

distance  between  atoms  i  and  j)  is  in  angstroms.  D  is  an  effective  dielectric 
constant  taken  to  be  2  in  all  cases. The  second  and  third  terms  in  eq.  (l) 
constitute  the  repulsive  and  attractive  van  der  Waals  interaction,  respectively, 
and  k  and  1  refer  to  different  types  of  atoms.  The  constants  employed  in  this 
calculation  are  given  in  Table  I  along  with  aromatic-aromatic  atomic  interaction 
parameters  for  completeness. 

A  hydrogen  bonding  interaction  term  utilizing  a  10-12  potential  replaces 
the  6-12  term  in  eq(l)  whenever  a  hydrogen  bonding  interaction  is  encountered. 
Specifically, 

'h.  ■  jj.  (%.  - 


in  which  the  sum  over  i  is  taken  only  over  the  hydrogen  bonding  atoms  H(N), 

H(0),  0,  or  N  in  the  first  molecule  and  the  sum  over  j  is  taken  over  the  other 
partner  atom  in  the  hydrogen  bond  belonging  to  the  second  molecule.  A'  and  B 
constants  depend  on  the  type  of  hydrogen  bond  formed,  in  the  sense  that  they 
depend  on  the  nature  of  the  three  atoms  involved. In  subsequent  reports 
we  will  use  eq.  (2)  in  the  calculation  of  the  binding  energy  and  geometry  of  the 
pyrazine  dimer,  pyrazine  ammonia  clusters,  and  pyrazine  water  clusters. 

For  comparison  purposes,  calculations  are  reported  for  both  LJ  and  exp-6 
potential  forms  on  benzene-methane  and  benzene-ethane  clusters.  The  calculated 


binding  energies  for  benzene-methane  using  these  two  potentials  are  540  cm 
and  589  cm~^ ,  respectively.  The  binding  energy  is  smaller  for  the  LJ  poten¬ 
tial  than  for  the  exp-6  potential  by  -  50  cm~^  but  the  geometries  are  virtually 

/  0  \ 

identical  (_+  0.05  Aj.  Similar  findings  obtain  for  the  benzene-ethane 
clusters:  two  geometries  are  predicted  using  either  potential^  with  binding 
energies  for  the  LJ  potential  about  50  cm  ^  below  those  for  the  exp-6 
potential.^  Fig.  1  displays  the  computer  drawings  of  these  configurations  for 
the  different  potentials.  These  two  potentials  give  indistinguishable  results 
for  binding  energies  and  geometry  estimates  to  within  the  experimental  accuracy 
In  all  other  respects  the  calculations  are  carried  out  as  previously  reported. 
III.  Results 

2-color  TOFMS  and  FE  data  are  reported  below  for  pyrazine  and  pyrazine 

clustered  with  n  =  1 ,  2,  3*  Clusters  containing  up  to  three 

solvent  molecules  are  identified  and  characterized.  For  each  cluster  both  0° 

0 

112 

and  higher  vibronic  transitions  observed:  experi¬ 

mental  data  and  potential  energy  calculations  will  be  presented  for  the  various 
clusters  formed. 

A.  Pyrazine 

*  1  1 
The  singlet  n->  ir  transition  of  pyrazine  <->  appears 

at  30,875.5  cm”\  The  10a\  6a\  and  10a^  states  appear  at  383,  583.7 

and  823  cm~^  above  the  0°  transition,  respectively.  The  pyrazine  spectrum 

(fig.  2)  is  identical  in  FE  or  2-color  TOFMS  and  is  the  same  as  has  been 

previously  published.  The  pyrazine  structure  used  in  the  calculations  is 

obtained  from  ref.  11. 

B.  Fyrazine-methane 


The  origin  of  the  0^H^N2(CH^)^  cluster  is  shifted  -33  cm 
from  the  pyrazine  0°  as  shown  in  Table  II  and  fig  3-  No  vdW  modes  are 


observed  built  on  this  origin.  2-color  TOFMS  signals  are  observed  for 

11  2 
the  10a  and  6a  transition  but  not  for  the  lOa"  transition  of 
0  0  o 

.  Thus,  the  binding  energy  for  the  cluster  is  between  684 
and  823  cm”^ . 

Based  on  comparison  with  other  aromatic-hydrocarbon  clusters^ and 
on  potential  energy  calculations,  the  minimum  cluster  energy  is  associated 
with  the  methane  molecule  coordinated  to  the  tt -system  of  the  aromatic  ring. 
Pig.  4  gives  a  top  view  of  the  CH^  above  the  pyrazine  plane.  The  calculated 
binding  energy  is  512  cm.  ^  We  expect  the  LJ  potential  to  give  somewhat 
lower  than  experimental  binding  energies,  as  presented  above  for  the  benzene 
test  case. 

The  C^H^N2(CH^)2  0°  spectrum  displays  two  peaks,  one  with  a  -67*9 
cm”^  shift  and  the  other  with  a  -39*1  cm  ^  shift.  These  data  are 
presented  in  fig.  3  and  Table  II.  The  larger  of  the  two  shifts  is  nearly  twice 
that  observed  in  the  C^H^N2(CH^)^  cluster:  the  implication  is  that  the 
second  methane  is  added  symmetrically  to  the  other  side  of  the 
^4^4^2^^^4^1  plane.  The  -39»1  cm  ^  shift  must  then  be 

associated  with  the  (anisotropic)  cluster  with  both  methanes  on  the  same  side 
of  the  aromatic  ir  -system. 

The  LJ  calculation  of  the  binding  energies  and  stable  geometries  for  the 
C^H^N2(CH^)2  cluster  gives  results  consistent  with  the  above  findings. 

One  calculated  cluster  has  methanes  on  either  side  of  the  pyrazine  plane  in 
mirror  image  positions  with  a  binding  energy  of  1023  cm  The  second 
*^4^4^2^^^4  ^2  features  both  methanes  on  the  same  side  of  the 

ring,  one  shifted  toward  a  nitrogen  atom  and  the  other  away  from  the  ring; 
this  configuration  gives  a  binding  energy  of  852  cm  ^ .  Fig.  4  shows  these 
geometries. 


The  C^H^N2(CH^)^  clusters  have  more  diffuse  spectra  as  is  typical 
for  larger  clusters.  The  spectra  suggest  two  configurations,  one  with  three 
methanes  on  one  side  of  the  aromatic  plane  (-35  cm”^ )  and  one  with  two 
methanes  on  one  side  and  one  on  the  other  (-71  cm”^).  Fig.  3  and  Table  II 
give  these  data. 

Calculations  were  not  performed  on  the  C^H^N2(CH^)^  clusters. 

C.  Pyrazine-Ethane 

The  0°  spectrum  of  ^4H4N2(C2Hg shows  three  features  of  roughly 

equal  intensity  which  can  best  be  attributed  to  three  different  geometries  (fig. 

5).  Two  of  these  features  are  separated  by  4  cm~^  and  appear  at  red  shifts  of 

56.3  and  52.5  cm~^ .  The  third  feature  has  a  shift  of  -41.1  cm"^ .  These 

are  listed  in  Table  III.  Some  cluster  vibronic  transitions  are  also 

identified  (tentatively)  to  the  high  energy  side  of  the  highest  energy 

C^H^N2(C2Hg)^  origin.  Vibronic  transitions  are  also  studied  for  this 

cluster  and  the  binding  energy  can  again  be  bracketed  as  reported  above  for 

1  2 

*^4^4^2^^^4^1  solvent  molecule. 

The  three  cluster  configurations  observed  in  the  spectrum  of 
C^H^N2(C2Hg^) ^  are  also  found  in  the  calculations.  These  are  presented  in 
fig.  6.  The  rotational  barrier  between  geometry  II  and  III  is  13  cm  ^  which 
is  probably  well  above  the  zero  point  motion  in  the  ethane  rotational  coordinate 
carrying  geometry  II  to  III  (fig.  6)  by  a  simple  clockwise  rotation.  Configura¬ 
tion  I  in  fig.  6  displays  the  ethane  molecule  with  its  long  axis  nearly  coin¬ 
cident  with  the  Z-axis  of  pyrazine  (X=0.00,  Y=  +  0.05  X).  We  have  assigned  this 
geometry  to  the  largest  of  the  three  cluster  shifts  because  ethane  interacts 
with  the  7T -system  most  strongly  in  this  geometry.  This  assignment  is  consistent 
with  the  benzene-ethane  geometry  with  the  larger  cluster  shift. ^  The  other  two 
geometries  find  the  ethane  in  asynmetric  positions  above  the  ring  with  the  long 


axis  of  the  ethane  more  or  less  parallel  to  the  plane  of  the  ring.  Configura- 

__  0 
tion  II  has  one  carbon  in  ethane  situated  5.84  A  above  the  C-N  bond  with  the 

other  carbon  5-59  X  above  the  ring  and  pointing  towards  the  diagonally  opposite 

C-N  bond.  Based  on  the  overall  configuration,  it  would  appear  as  though  some 

"hydrogen  bonding"  (.i.e.,  hydrogen-nitrogen  preferred  orientation)  is  important 

here,  even  though  a  specific  hydrogen  bonding  potential  is  not  employed  in  the 

calculation.  In  the  third  geometry,  configuration  III  in  fig.  6,  the  ethane  is 

also  in  an  asymmetric  position:  here  one  carbon  is  3.81  K  above  C-C  bond  and 

the  other  carbon  is  3.56  X  above  the  ring  pointing  towards  the  other  C-C  bona. 

Once  again  the  hydrogens  are  oriented  towards  the  nitrogen  atoms. 

The  spectrum  of  pyrazine  with  two  ethanes  shows  three  features  which  can 

be  associated  with  C  .H  .N^(C,,H.  )„  0°  transitions.  Two  of  these 

4  4  2  2  6  2  o 

transitions  are  intense  and  appear  at  -65.0  and  -57.6  cm"^  from  tr.e  pyrazine 

0°  transition.  The  third  feature  is  identified  at  -41 .8  cm”^  from  the 
0 

pyrazine  0®  and  is  quite  weak.  The  configurations  to  which  we  have  assigned 
these  features  are  presented  in  fig.  7.  Note  that  none  of  these  features 
evidences  an  additive  shift  with  respect  to  the  ^4H4N2(C2H^)^  spectra 
and  thus  the  isotropic  or  symmetric  2^^ 2^6^ 2  configurations  are  not 

responsible  for  the  major  observed  intensity.  It  is,  however,  possible  that 
the  very  weak  and  broad  feature  at  -78  cm~^  may  be  associated  with  an  additive 
shift  cluster. 

Potential  energy  calculations  give  four  possible  low  energy  asymmetric 
or  anisotropic  geometries  for  the  C^H^N^Cc^H^)^  clusters,  as  depicted  in 
fig.  7.  In  each  configuration  the  ethanes  are  on  the  same  side  of  the  aromatic 
ring,  as  mandated  by  the  absence  of  additive  cluster  shifts.  Two  distinctive 
orientations  of  the  ethanes  above  the  ring  are  noted.  In  the  first  case, 
configuration  I  and  II,  one  ethane  is  perpendicular  to  and  near  the  center  of 
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the  aromatic  (XY)  plane  while  the  other  ethane  is  parallel  to  the  plane  and 
displaced  from  the  ring  center.  These  two  geometries  are  rather  similar  in 
binding  energy  and  probably  would  have  similar  cluster  shifts.  The  second  case 
(III  and  IV)  has  both  ethanes  parallel  to  the  ring  in  various  positions.  The 
binding  energies  associated  with  configurations  III  and  IV  (fig.  7)  are  also 
quite  similar.  Since  specific  assignments  of  spectra  to  geometry  can  only  be 
made  in  this  instance  through  comparison  with  benzene-ethane  clusters,  we  post¬ 
pone  the  discussion  of  such  a  correspondence  until  the  next  section. 

Since  the  major  spectra  for  clusters  are  founa  at  roughly 

the  same  shifts  as  the  spectra  for  ^4H^N2(C2Hg)^  2  clusters,  we  again  conclude 

that  the  three  C„H^  molecules  are  all  on  the  same  side  of  the  ring.  Note  that 
td  o 

the  base  line  has  now  built  up  under  the  main  features,  probably  due  to  the  large 
number  of  possible  configurations  for  the  ethanes  not  directly  coordinated  to  the 
pyrazine  ring. 

D.  Pyrazine-Propane 

The  pyrazine-propane  system  shows  results  not  very  different  from  those 
presented  above  and  similar  to  those  found  for  benzene-propane.^  The 
spectrum  of  C^H^N2(C^Hg) ^  gives  only  one  peak  with  a  red  shift  with 
respect  to  the  pyrazine  origin  of  60.3  cm  ^  (see  fig.  8  and  Table  IV).  No 
other  features  are  observed  for  this  cluster.  The  calculations  associate  the 
spectrum  with  a  geometry  in  which  the  propane  lies  horizontally  on  the  pyrazine 
ring  and  roughly  3-7  X  above  it  as  is  depicted  in  fig.  9-  This  configuration 
has  a  calculated  binding  energy  of  883  cm  ^ . 

For  pyrazine  with  two  propanes  we  again  find  no  additive  shifts  and  thus 
conclude  both  propane  molecules  are  on  the  same  side  of  the  ring.  While  the 
spectrum  is  in  general  weak,  a  feature  is  evident  at  -10if.3  cm  ^  and  two 
others  are  clearly  observed  with  smaller  shifts  (see  fig.  8  and  Table  IV). 
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A  number  of  potential  energy  minima  exist  for  the  C^H^N2(CjHg)2  system: 
we  have  not  explored  the  entire  surface  for  this  cluster.  One  configuration 
with  the  propanes  arranged  as  shown  in  fig.  9  is  stable  with  a  1012  cm”^  bind¬ 
ing  energy.  Many  other  configurations  are  undoubtedly  stable  for  this  cluster. 
Table  V  gives  a  list  of  the  calculated  and  experimental  binding  energies. 

IV.  Discussion 

The  combination  of  2-color  TOFMS  (mass  and  optical  selectivity)  and 
potential  energy  calculations  can  be  employed  to  determine  cluster  geometry 
and  nucleation  processes.  A  comparison  between  pyrazine-hydrocarbon  and 
benzene-hydrocarbon  clusters  will  be  presented  in  this  section  along  with  a 
brief  presentation  of  the  nucleation  processes  in  the  beam.  Because  of  the 
overall  similarity  between  the  calculations  and  data  reported  in  this  work 
and  those  found  in  ref.  1-6  for  the  benzene,  toluene  and  aniline  clusters, 
we  present  a  short  summary  of  the  major  points  made  in  these  published  papers. 
The  conclusions  are  as  follows: 

1 .  spectral  cluster  shifts  are  intimately  connected  with  the  geometry  of 
the  cluster; 

2.  the  size  of  the  shift  is  governed  by  the  solvent  molecule  overlap  with 
the  aromatic  tt- system  of  the  solute; 

3*  the  size  of  the  shift  scales  with  solvent  polarizability  and  the  change 
in  solute  polarizability  upon  excitation  -  hence,  the  shift  is  also 
dependent  upon  solvent  orientation  with  respect  to  the  solute; 

4.  additivity  of  cluster  shifts  for  equivalent  structural  sites  on  the 
solute  molecule  is  maintained  to  within  ^2  cm~^  for  benzene,  toluene 
and  aniline  solvated  by  methane,  ethane  and  propane; 


5.  calculations  of  binding  energies  with  the  exp-6  atom-atom  potential  gives 
binding  energies  to  within  the  experimental  accuracy; 


6.  according  to  the  exp-6,  L-J  benzene  comparison  calculation,  the  calculated 
LJ  binding  energy  is  expected  to  be  ca.  50  cm~^  below  the  calculated 
exp-6  binding  energy; 

7-  both  potentials  give  geometries  that  are,  to  within  calculations! 
accuracy,  identical; 

8.  the  calculated  geometries  agree  with  tne  benzene  experimental  geometries 
which  are  uniquely  identified  based  on  symmetry  consideration;  and 
9-  in  all  instances  for  which  comparisons  are  possible,  the  calculations 
agree  with  the  experimental  findings  and  vice  versa. 

A.  Pyrazine  -  Methane 

« 

Pyrazine  solvated  by  methane  produces  a  set  of  n  ->•  it  spectra 
that  are  in  general  simple  to  interpret.  The  C^H^N2(CH^)^  complex 
evidences  only  one  transition  for  each  vibronic  feature  of  pyrazine  observed. 

No  vdW  mode  vibronic  additions  to  the  cluster  0°,  10a\  6a^  or  10a^  are 
observed.  We  attribute  the  absence  of  Av  =  ^1  vibronic  transitions  to  the 
similarity  of  the  ground  and  excited  state  vdW  potential  energy  surfaces. 

These  surfaces  are  more  nearly  identical  than  those  of  any  other  system  we 

have  studied. This  similarity  can  be  rationalized  based  on  the  nature 

* 

of  the  electronic  transition.  In  a  w  -►  tt  transition,  an  electron 
involved  in  the  ground  state  binding  of,  for  example,  benzene  and  methane  is 

removed  from  a  bonding  orbital  and  placed  in  an  antibonding  orbital.  However, 

* 

the  n  -+  IT  transition  removes  a  nonbonding  electron,  probably  not  very  much 
involved  in  the  overall  pyrazine  methane  interaction,  and  places  it  in  an  anti¬ 
bonding  orbital  again  not  heavily  involved  in  the  cluster  bonding.  Thus 
the  ground  and  excited  state  potential  energy  surfaces  must  be  nearly 
superimposable ,  even  though  the  excited  state  well  is  -40  cm  ^  deeper.  More¬ 
over,  this  argument  suggests  that  the  and 

clusters  should  have  similar  shifts  because  both  ground  and  excited  state 
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IT -systems  are  nearly  the  same.  The  methane  in  a  pyrazine-methane  complex 
coordinates  with  the  aromatic  it -system  in  much  the  same  manner  as  found  for 

V6<™4^- 

The  C^H^N2(CH^)2  system  is  again  very  similar  to  the  CgH^(CH^)2 
system.^  Both  isotropic  and  anisotropic  configurations  are  observed.  The 
intensity  ratio  of  the  isotropic  to  anisotropic  clusters  is  roughly  2:1,  again 
as  found  for  the  benzene  system.  The  second,  asymmetric  methane  does  not  over' 
lap  the  ring  and  thus  contributes  little  to  the  cluster  shift  which  is  nearly 
the  same  as  found  for  C^H^N2(CH^)^ .  The  binding  energy  for  the  symmetric 
cluster  is  1023  cm~^  and  for  the  asymmetric  cluster  is  852  cra”^ .  , 

on  the  other  hand,  has  both  symmetric  and  asymetric  clusters  with  nearly  the 
same  binding  energy  (approximately  1180  cm  ^). 

B.  Pyrazine-Ethane 

The  cluster  can  have  three  distinct  geometrical 

arrangements  with  roughly  equal  probability  and  similar  binding  energies  (fig. 
6),  whereas  only  two  configurations.  This  difference 

must  be  a  direct  consequence  of  the  presence  of  the  nitrogen  atoms  in  the  aro¬ 
matic  ring.  The  cluster  with  the  largest  spectral  shift  (-56*5  cm  has  been 
assigned  as  the  one  for  which  the  ethane  is  perpendicular  to  the  pyrazine  ring. 
In  C,H, (C„H^),  this  geometry  has  a  shift  of  -57.2  cm  ^ .  In  both  instances 
this  configuration  has  the  smallest  binding  energy.  Weak  vibronic  Av  =  j^1 
transitions  are  observed  for  the  C,H,N„(CoH£ ),  clusters  and  thus  the  vdW 
potentials  for  the  ground  and  excited  states  are  somewhat  different. 

We  can  assign  the  two  parallel  orientation  clusters  to  specific  shifts 
(-52.5  and  -41.1  cm”^ )  in  the  same  manner  suggested  for  the  aromatic  hydro¬ 
carbons:  the  cluster  with  the  larger  overlap  between  the  solute  tt -system  and 
the  solvent  has  the  larger  cluster  shift.  Thus  the  geometry  (II  in  fig.  6) 


with  the  C_H^  oriented  toward  the  C-N  bond  is  associated  with  the  -52.5  cm 
d  b 

shift  and  the  geometry  with  the  oriented  toward  the  C-C  bond  is  associated 

with  the  -41.1  cm”^  shift.  This  correlation  is  at  least  consistent  with  all 
other  systems  we  have  studied  to  date. 

Three  distinct  features  are  observed  for  the  C .H )„ 

0,  o,  d  d  b  d 

clusters.  Figs.  5  and  7  present  the  spectra  and  calculated  geometries. 
Calculated  configurations  I  and  II  are  quite  similar,  as  are  calculated  con¬ 
figurations  III  and  IV.  In  either  grouping  only  one  ethane  is  well  coordinated 
with  the  aromatic  ring;  in  the  first  set  the  coordinated  ethane  is  perpendic¬ 
ular  to  the  ring  and  in  the  second  it  is  parallel.  The  nitrogen  atoms  seem  to 
play  a  central  role  in  these  geometries  because  for  2 

ethanes  are  always  perpendicular  to  the  ring. 

We  suggest,  baaed  on  arguments  and  correlations  made  previously,^  that 
configuration  I  and  II  can  be  associated  with  spectral  cluster  shifts  of  -57.6 
and  -65*0  cm”^ .  Configuration  II  can  be  assigned  to  the  -65*Ocm  ^  feature 
for  two  reasons:  the  perpendicular  ethane  has  slightly  more  overlap  with  the 
ring  in  this  geometry  and  this  geometry  has  the  slightly  lower  binding  energy. 
The  shift  at  -41 .1  cm~^  is  attributed  to  the  configuration  IV  since  one  ethane 
has  an  overlap  with  the  ring  similar  to  the  C^H^N2(C^H^)^ case.  Configur¬ 
ation  III  is  not  readily  assigned  to  any  of  the  observed  features. 

C.  Pyrazine-Propane 

The  geometry  in  this  instance  for  C^H^NgCC^Hy) ^  is  very  similar 
to  that  found  for  C^H^(C^Hg)^.  Apparently  the  nitrogen  atoms  do  not 
affect  the  orientation  of  the  propane  above  the  pyrazine. 

D.  Nucleation 

Both  symmetric  (isotropic)  and  asymmetric  (anisotropic)  clusters  are 
possible  for  the  solute(solvent)  systems.  The  isotropic  clusters  have  a 


advent  molecule  on  either  side  of  the  ring  and  the  anisotropic  clusters  have 
both  solvent  molecules  on  the  same  side  of  the  ring.  Both  types  of  solute- 
(solvent)^  clusters  can  be  formed  through  homogeneous  nucleation,  ’  the 
addition  of  one  solvent  molecule  at  a  time  to  the  cluster.  Inhomogeneous 
nucleation,  the  addition  of  two  molecules  to  the  cluster  simultaneously,  can 
only  produce  anisotropic  clusters.  We  have  previously  shown  that  the 
supersonic  molecular  jet  conditions  (i.e.,  backing  pressure,  nozzle  diameter, 
pumping  speed,  concentration  ratios)  used  in  our  laboratory  deplete  the  beam  of 
solvent  monomer  species. 

The  C^H^N2(CH^)2  cluster  spectrum  show  isotropic  and  anisotropic 
clusters  in  roughly  similar  amounts.  Pyrazine  expanded  with  ethane  or 
propane,  on  the  other  hand,  gives  only  anisotropic  clusters.  These  concen¬ 
tration  ratios  are  a  consequence  of  the  nucleation  processes  in  the  beam  as  we 
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have  demonstrated  for  benzene  and  toluene.  ’  A  complete  discussion  of  these 
mechanisms  can  be  found  in  our  previous  papers. 

V.  Conclusion 

Pyrazine-hydrocarbon  clusters  are  in  general  similar  to  benzene 
hydrocarbon  clusters:  the  presence  of  nitrogen  atoms  in  the  aromatic  ring  has 
a  relatively  small  although  discernable  effect  on  the  overall  solute-solvent 
interaction.  In  particular,  solvent  hydrogen  atoms  preferentially  orient  toward 
the  solute  nitrogen  atoms.  While  some  minor  differences  do  exist,  the  cluster 
spectral  shifts  for  the  two  series  are  quite  similar  as  well. 

Calculations  of  geometries  and  binding  energies  using  either  the  LJ  or 
an  exp-6  potential  produce  nearly  identical  results  for  benzene  and  pyrazine 
species.  The  results  based  on  either  potential  form  are  in  agreement  with 
experiments  and  are  self  consistent  for  all  species  thus  far  studied. 

Inhomogeneous  nucleation  predominates  as  the  main  process  in  the  formation 
of  solute(solvent)_  clusters. 
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Table  I 


Parameters  for  the  energy  expression  in  the  computer  modeling 


E..=1.16xl(?  ‘’iV2r  + 
ij  ij  ij 


/  -1  .ol2\  /  -1  .o6\ 

I  cm  A  I  /cm  A  | 

\  mole  /  y  mole  / 


-  ^Wr^ 

ij 


A1 ip ha  t ic-Al ipha  t ic 


c-c 

3.168 

X 

10® 

1.295 

X 

10® 

C-H 

4.182 

X 

10^ 

4.397 

X 

lo"* 

H-H 

4.926 

X 

10® 

1.591 

X 

10^ 

■Aromatic 

C-C 

2.711 

X 

10® 

1.517 

X 

10® 

C-H 

4.217 

X 

10^ 

4.397 

X 

10^ 

H-H 

4.979 

X 

10® 

1.590 

X 

lo"* 

C-N 

2.077 

X 

10® 

1.348 

X 

10^ 

H-N 

2.500 

X 

10^ 

4.527 

X 

10^ 

iromatic 

C-C 

2.285 

X 

10® 

1.781 

X 

10® 

C-H 

3.462 

X 

10^ 

5.217 

X 

10^^ 

H-H 

5.028 

X 

10® 

1.591 

X 

10^ 

C-N 

1.728 

X 

10® 

1.575 

X 

10® 

H-N 

2.523 

X 

10^ 

4.527 

X 

10^ 

N-N 

1.312 

X 

10® 

1.403 

X 

10® 

Aromatic-Aliphatic 


Table  II 


Observed  peaks  in  the  spectra  of  pyrazine-methane  clusters. 


Species 

Energy_j 
(vac .  cm  ) 

Energy  Relative 
to  Pyrazlne  oS 
(cm-1) 

Energy  Relative 

to  Cluster  0° 
(cm-1)  o 

Assignment 

30842.9 

-33.1 

0 

o  o 
o 

C4H^N2^^”4^2 

30808.1 

-67.9 

0 

iso  0*^ 
o 

30836.9 

-39.1 

0 

aniso  0^ 
o 

C4H4N2(CH^)3 

30805.0 

o 

• 

1 

b 

o  o 
o 

30841.0 

-35.0 

0 

o  c 
o 

-1 


a)  The  ~67.9  cm  shift  is  associated  with  the  two  methanes  added  symmetrically 
above  and  below  the  pryazine  ring  (Isotropic),  as  shown  in  fig.  4.  The  -39.1 
cm”^  shift  is  associated  with  the  anisotropic  configuration,  as  shown  in  fig.  4. 


Table  III 


Observed  peaks  In  the  spectra  of  pyrazine-ethane  clusters 


Species 


Energy 
(vac.  cm  ) 


Energy  Relative 

to  Pyrazine  0° 
(cm-1)  o 


Energy  Relative 

to  Cluster  Oq 
fcm-1) 


Assignment 


30819.7 

-56.3 

0 

I  0 

o 

30823.5 

-52.5 

0 

II  0° 

30834.9 

-41.1 

0 

III  0° 

30839.5 

-36.5 

16.0 

II 

O 

30845.3 

-30.7 

10.4 

III  A^ 
o 

30854.3 

-21.7 

30857.9 

-18.1 

W2<-W2 

30811.0 

-65.0 

0 

II  0*" 
o 

30818.4 

-57.6 

0 

I  0° 

30834.2 

-41.8 

0 

III  0° 

W2<V6’3 

30803.3 

-72.7 

30810.5 

-65.5 

30816.7 

-59.3 

30821.5 

-54.5 

30829.9 

-46.1 

30835.2 

-40.8 

a).  V  »  a  vdW  stretch 

A  =  a  vdW  bend 

I,  II,  III  =  represents  different  configurations 


Assignment 


Energy  Relative  Energy  Relative 

Energy_^  to  Pyrazine  to  Cluster 

Species  (vac.  cm  ■*■)  Og  (cm“^)  0°  (cm“^) 

s 


I*'  *' 


SV2<S«8>1 

30815.7 

-60.3 

0 

0° 

Wl(S«8>2 

30773.7 

-102.3 

0 

o  o 
o 

30807.5 

-68.5 

0 

o  o 
o 

30826.5 

-49.5 

0 

0° 

o 


^  a)  Different  possible  configurations  of  C.H.N. (C_H„)„  give  rise  to  the  shifts  at  -102,3, 
Ij  -68.5,  and  -49.5  cmT^  See  fig.  9  for  an  example. 


Table  V 


Binding  energies  for  C,H, 

O  D 

and  C.H.N^ 

4  4  2 

with  various  hydrocarbons. 

Species 

Configuration 

Calculated  (cm  Experimental 

LJ^  EXP-6  (cm”^) 

CeHeCCH^)! 

-540 

-589 

I 

-598 

-654 

S«6<S«6>1 

II 

-705 

-778 

C4H4N2(CH^)j 

-512 

584  <  B.E.  <  823 

iso 

-1023 

C,H^N2(CH^)2 

aniso 

-852 

W2<S«6h 

I 

-564 

584  <  B.E.  <  323 

V4N2<S«6h 

II 

-668 

W2<SVl 

III 

-672 

V4N2^S»6^2 

I 

-1089 

'^4V2<S»6>2 

II 

-1086 

S»4^2^S»6^2 

III 

-1198 

W2^''2»6>2 

IV 

-1194 

W2<S«8h 

-883 

W2(S«8h 

-1012 

a)  The  LJ  potential  gives  lower  binding  energies  than  the  exp-6  potential,  although  the 
geometrical  configurations  calculated  with  both  potentials  are  Identical. 

Figure  Captions 


Comparison  of  the  geometries  and  binding  energies  of  C^H^(CH^)  and 

clusters,  based  on  calculation  using  Lennard-Jones  (LJ)  and 

DO  o 

exponential-6  (exp-6)  potentials.  Geometries  are  shown  looking  down  the 
z-axis.  The  coordinate  axes  are  displayed  at  the  top  of  the  figure. 


Figure  2 

Time  of  flight  mass  spectra  (TOFMS)  of  pyrazine,  with  the  energy  scale 

relative  to  the  0°  transition  (30875.5  cm  ^ ) .  The  three  main  peaks 

2  -1 

to  higher  energy  are  10a,  6a,  and  10a  at  383,  584,  and  823  cm,  respectively. 


Figure  3 

TOFMS  of  C.H,N„(CHj,,  C.H.N-(CHJ„,  and  C.H.N»  (CHj,  in  the 
4  4  2  4  1  44  2  4  2  442  43 

region  of  the  pyrazine  origin.  The  uppermost  trace  is  the  optogalvanic  signal 
used  for  calibration  purposes:  on  this  and  the  other  spectra  presented,  the 
feature  appears  at  6508.3259  X  (vacuum).  The  pyrazine  origin- falls  at  0  cm  ^ 
on  the  scale. 


Minimum  energy  configurations  and  binding  energies  for  C^H^N^  (CH^)^  and 
^4^4^2^^^4^2  obtained  using  a  LJ  potential  calculation. 


Figure  5 

TOFMS  of  C^H^N2(C2H^)^ ,  C^H^N^(C2H^)^,  and  C^H^N2(C2Hg)^  in 

the  region  of  the  pyrazine  origin  (displayed  as  0  on  the  energy  scale)- 


Minimum  energy  configurations  and  binding  energies  for  as 

obtained  using  a  LJ  potential  calculation. 

Figure  7 

Minimum  energy  configurations  and  binding  energies  for  as 

obtained  using  a  LJ  potential  calculation. 

Figure  8 

TOFMS  of  and  in  the  region  of  pyrazine 

origin. 

Figure  9 


Minimum  energy  configurations  and  binding  energies  for  C^H^N2(CjHg)^  and 

using  a  LJ  potential  calculation. 
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